Magnetic Phase Control in Monolayer Films by Substrate Tuning 
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We propose to tailor exchange interactions in magnetic monolayer films by tuning the adjacent 
non-magnetic substrate. As an example, we demonstrate a ferromagnetic-antiferromagnetic phase 
transition for one monolayer Fe on a Ta^Wi-^OOl) surface as a function of the Ta concentration. 
At the critical Ta concentration, the nearest-neighbor exchange interaction is small and the magnetic 
phase space is dramatically broadened. Complex magnetic order such as spin-spirals, multiple-Q, 
or even disordered local moment states can occur, offering the possibility to store information in 
terms of ferromagnetic dots in an otherwise zero-magnetization state matrix. 

PACS numbers: 75.70.Ak, 71.15.Mb 



Magnetic systems play a central role in today's infor- 
mation technology and our ability to control and tailor 
their properties may open new vistas to future device 
concepts. Materials structured on a nanometer scale such 
as atomically thin films proved to be a rich field for novel 
magnetic properties. While our understanding of mag- 
netic systems has tremendously increased over the past 
25 years, to control magnetic order in a specific system 
and tailor materials with desired magnetic properties re- 
mains the grand challenge of research in magnetism. 

So far the attempts to tune the magnetic state of sur- 
faces and ultra-thin films have focussed on alloy for- 
mation where the concentration of the magnetic com- 
ponents is altered to optimize the magnetic proper- 
ties [l|, d, 0, ID, Si- In this letter, we propose a com- 
pletely different route. Based on the surprising observa- 
tion of an antiferromagnetic (AFM) order of one mono- 
layer (ML) Fe on W(001) and the prediction of the ferro- 
magnetic (FM) order on Ta(OOl) |7J, we suggest to tunc 
magnetic interactions in ultra-thin films by modifying 
only the band-filling of the substrate, via the formation 
of a Ta-W alloy, without altering or diluting the magnetic 
monolayer itself. 

We employ first-principles calculations to show that 
the nearest-neighbor exchange interaction in one ML Fe 
on the (001) surface of a Ta^Wi-j alloy can be con- 
tinuously tuned from FM to AFM coupling by varying 
the Ta concentration x. At the substrate composition 
of small nearest-neighbor exchange interaction, we find 
that higher order spin interactions beyond the Heisenberg 
model, such as biquadratic or four-spin interactions, may 
stabilize complex non-collinear magnetic structures. In 
this case, we also consider the role of chemical disorder 
which might prevent any stable magnetic order due to 
the small energy scale involved and lead for example to a 
spin-glass. The substrate turns out to be a tuner also for 
a magnetic order-disorder phase transition. At the corre- 
sponding Ta concentration, a highly frustrated material 



is formed and can be used to store information in the 
form of FM dots in a zero-magnetization state matrix, 
opening the way to a new class of material for magnetic 
storage devices. 

While we consider a single model system, which allows 
to isolate magnetic effects from structural and chemi- 
cal ones, depending on surface orientation, substrate el- 
ement, and overlayer, a rich variety of similar systems 
are possible. For example, metallic magnets with small 
exchange coupling have been recently reported [3, [3] . 

We have determined the electronic and magnetic 
properties of one ML Fe on the (001) surface of 
Ta^Wi-a; by performing density-functional theory calcu- 
lations in the generalized-gradient approximation to the 
exchange-correlation functional [13]. The substitutional 
Ta^Wi-a; random alloy has been modelled in the spirit of 
the virtual crystal approximation (VCA) 11] by a sub- 
strate of fictitious atoms with fractional atomic numbers 
related to the Ta composition x, ranging linearly between 
73 (Ta) and 74 (W). The corresponding fractional elec- 
tronic charge preserves charge neutrality and accounts for 
the variation of the band-filling originating from alloying. 
Vegard's law was adopted to interpolate between the lat- 
tice constants of Ta (3.301 A) and W (3.165 A). A fixed 
surface relaxation of 18% was assumed, corresponding to 
the relaxation of the FM Fe monolayer on pure W(001). 
Based on additional studies on the effect of the lattice 
constant on the magnetic properties, the tiny deviations 
from Vegard's law llj and the deviation of the relaxation 
from the average value [13] can be safely neglected. The 
calculations have been carried out with the full-potential 
linearized augmented plane wave (FLAPW) method in 
film geometry, as implemented in the fleur code [1 31 ) - 
Spin spirals have been computed in the p(l x 1) unit cell, 
exploiting the generalized Bloch theorem [14j | . The com- 
putational parameter were chosen according to Ref. [7j. 

W and Ta are adjacent elements of the periodic ta- 
ble with similar properties. Both crystallize in the bcc 
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FIG. 1: (color online) (a) Unit cell sketch of the investigated 
magnetic configurations, (b) Total energy, (c) magnetic mo- 
ment, and (d) exchange constants for f ML Fe/Ta^Wi-a; as 
calculated by the FLAPW method. 



structure with comparable lattice constants. W has one 
d electron more than Ta. While one ML Fe exhibits a 
c(2 x 2)-AFM state, Fig. QJa), if grown on W(001), it 
is FM on Ta(001) [7J. Hence, the magnetic configura- 
tion in the Fe layer is related to the substrate d-band 
filling that affects the position of the substrate d-band 
relative to the Fermi energy (Ep) and controls by hy- 
bridization the position of the Fe states at Ep. For 3d 
metals with large magnetic moments it is well-established 
that in good approximation the magnetic configuration 
with the lowest minority density of states at Ep exhibits 
the lowest energy (e.g. compare Fig. 3 and 4 in Ref. Ila). 
By forming TaaWi—x alloys with different compositions 
we can adjust the 5d-band filling or the position of the 
minority Fe d-state at Ep and continuously investigate 
the transition between the AFM and FM state of the 
Fe overlayer in Fe/W(001) and Fe/Ta(001), respectively. 
This is expected to be experimentally feasible due to the 
full miscibility of W and Ta for every concentration [n| . 

This scenario is presented in Fig. [TJ We initially focus 
on three collinear configurations, Fig.QJa): c(2x2)-AFM, 
p(2 x 1)-AFM, also referred to as row-wise AFM (RW- 
AFM), and FM. The magnetic ground-state is obtained 
from the total energy, Fig. []Jb), as a function of the 
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FIG. 2: (color online) Phase diagrams for 2D systems on a 
square lattice based on the Heisenberg model, (a) If Ji and J2 
dominate, only collinear states are possible, (b) If Ji is small 
compared to J3, the FM and the c(2 x 2)- AFM configurations 
are degenerate and spin spirals with q along T — M can occur. 



Ta concentration, x, relative to the FM state. At small 
Ta concentrations up to about 25%, the Fe monolayer 
exhibits a c(2 x 2)-AFM order, while at large x, roughly 
beyond 70%, the ground state is FM. In the intermediate 
range the RW-AFM state is energetically favorable ■ 

In order to analyze these ab-initio calculations in terms 
of exchange interactions, we can map the results onto a 
classical Heisenberg Hamiltonian, H = — J2i<j Jij s i ' s j> 
where Jij is the exchange interaction between spins at 
lattice sites i and j pointing in the direction of the unit 
vectors s.; and Sj, respectively. Such a model with fixed 
spin values is legitimated by the weak dependence of 
the magnetic moment on the Ta concentration and mag- 
netic order, Fig. [He). For a monolayer on a square lat- 
tice the nearest neighbor (nN) and next-nearest neighbor 
(nnN) exchange constants J\ and J2, respectively, can 
be extracted from the energies of the FM, c(2 x 2)-, and 
p(2 x 1)-AFM states. 

The nN exchange constant varies as a function of the 
Ta concentration, Fig. [Hd), and mimics the energy of 
the c(2 x 2)-AFM configuration, Fig. [T{b). In contrast, 
we observe only a weak dependence for the nnN coupling 
which is on the order of Ji ~ —5 meV. The strength 
of the typically dominating nN exchange interaction in 
a magnetic film can therefore be tuned by selecting the 
proper substrate composition, without requiring struc- 
tural or chemical modifications of the film. This opens 
the way to a dramatically broader magnetic phase space, 
as we will see in the following. 

The possible magnetic states in a two-dimensional (2D) 
film on a square lattice can be studied in the magnetic 
phase diagrams derived within the Heisenberg model, 
Fig. [2j For dominating nN exchange interaction only 
two magnetic states can occur, either the FM (Ji > 0) 
state or the c(2 x 2)-AFM one {J x < 0), Fig.^a). Even 
in the presence of nnN exchange the magnetic states are 
still limited to collinear solutions. However, if J\ is neg- 
ligible or at least comparable to interactions beyond sec- 
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FIG. 3: (color online) (a) Magnetic moment and (b) energy 
dispersion for spin spirals of an Fe UML at the W(001) lattice 
constant (dotted green line), 1 ML Fe/W(001) (solid red line), 
and 1 ML Fe/Ta. 42 W. 5 8 (001) (dashed blue line), (c) Disper- 
sion for 1 ML Fe/Ta.42W,58 (001) on a larger scale. Symbols 
denote ab-initio calculations and lines are Heisenberg fits up 
to eight neighbors. 



ond nearest neighbors, complex non-collinear magnetic 
phases appear, Fig. [2jb). For example, a considerable 
value of J3 leads to so-called spin-spiral states. Flat spin 
spirals are the general solution of the Heisenberg Hamil- 
tonian and are characterized by a wave vector q. For a 
given q, the magnetic moment of an atom at site R points 
in the direction s (R) = (cos (q • R) , sin (q • R) , 0). 

Based on these phase diagrams, the Ta concentration 
of very small J\, i.e. close to x = 42%, is particularly in- 
triguing, as interactions beyond nN will determine the 
magnetic ground state and non-collinear solutions are 
likely @, This suggests to extend the calculations 
by including spin spirals for vectors q along the high- 
symmetry lines of the irreducible 2D Brillouin zone. The 
high symmetry points represent the previously discussed 
collinear states: T, M, and X correspond to the FM, 
c(2 x 2)-AFM, andp(2 x 1)-AFM state, respectively. 

In order to emphasize the interplay of overlayer- 
substrate interaction and magnetic order, we study the 
energy dispersion E(q) for one ML Fe in different envi- 
ronments, Fig. 03 The unsupported monolayer (UML) 
of Fe(001) at the W lattice constant is an example of a 
system dominated by Ji, i.e. the FM state is the ground- 
state, the c(2 x 2)-AFM state is unstable, and the RW- 
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FIG. 4: Total energy of the configuration depicted in the inset 
for 1 ML Fe/Ta.42W.58 (001), as a function of the angle a. 



AFM state is metastable, with large energy differences 
among them. From a fit based on the Heisenberg model, 
we obtain the exchange constants and find J\ = 17 meV 
and J\ 3> Ji, for i > 1, in good agreement with the frozen 
magnon calculation by Sandratskii et al. [l8| ■ For 1 ML 
Fe/W(001) the hybridization with the substrate modi- 
fies the electronic structure. The dispersion is reversed, 
with the minimum at the M point. This implies that 
Ji is again the leading term, but with a negative sign: 
Ji = -26 meV. 

However, for 1 ML Fe/Ta.4 2 W.5 8 (001) the dispersion 
is strikingly different from the two previous cases. Since 
3\ ~ for this particular substrate composition, the 
FM and c(2 x 2)-AFM states are degenerate and the en- 
ergy scale is strongly reduced, Fig.OJb). Interestingly, a 
metastable spin-spiral state is found along the T-M direc- 
tion, 20 meV/Fe-atom lower than the FM one, Fig.[3jc). 
The RW-AFM state is the global energy minimum. The 
fit of the dispersion reveals that the leading interaction 
is J2 = —6 meV and the system is in the lower left part 
of the p(2 x 1) region in the phase diagram of Fig. [2fb). 

Due to the surface symmetry, there are two equivalent 
X states in the Brillouin zone (see inset of Fig. 02b)), 
corresponding to two degenerate p(2 x 1)-AFM configu- 
rations with perpendicular orientations of the ferromag- 
netically coupled rows. Any superposition of these two 
spin spirals, a so-called multi-Q state, is a degenerate 
solution of the Heisenberg Hamiltonian. However, the 
degeneracy with the p(2 x 1)-AFM state can be lifted by 
higher order interactions beyond the Heisenberg model, 
such as the four-spin and the biquadratic ones [191 ]. 

Such interactions are implicitly included in the ex- 
change correlation potential, and we can evaluate their 
magnitude from first principles. We performed calcu- 
lations in the p{2 x 2) unit cell rotating the moments 
on all sites by an angle a as depicted in the inset of 
Fig. | a = 0° corresponds to the RW-AFM state while 
for a = 45° we obtain the 2Q-state, a 2D non-collinear 
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FIG. 5: (color online) (a) Total energy and (b) magnetic mo- 
ment for 1 ML Fe/Ta^Wi-a; (001) as a function of the Ta 
concentration x, as obtained with the TB-LMTO method. 



structure with perpendicular adjacent moments. Since 
all these states are degenerate within the Heisenberg 
model, the total energy difference depends only on higher 
order interactions. We find that non-collinear states gain 
energy on the order of 5 meV due to these interactions, 
the minimum being at a = 31°. The fitting revealed 
that even terms beyond the biquadratic and four-spin 
interactions are present in this system. The moment ar- 
rangement can be slightly modified by the magnetocrys- 
talline anisotropy. E.g., an out-of-plane easy axis with 
an anisotropy energy on the order of 2 meV, similar to 
that of 1 ML Fe/W(001) 0, would decrease a from 31° 
to about 20°, based on Fig. [5J but would not prevent the 
non-collinear order. 

The small energy scale found in the spin spiral 
as well as in the 2Q-state calculation for 1 ML 
Fe/Ta.42W.58 (001) indicates that only a small amount 
of energy is required to rotate the magnetic moments 
and that there is a competition between several mag- 
netic interactions. This suggests a frustrated system. 
Under these circumstances, the disordered local moment 
(DLM) state with local moments pointing in random di- 
rections and zero net magnetization needs to be con- 
sidered. We have evaluated its energy using the tight- 
binding linear muffin-tin orbital (TB-LMTO) method in 
the atomic sphere approximation within the framework 
of the coherent-potential approximation (CPA) [2(| 21 1. 
Concerning the magnetically ordered states, we found 



that for small, intermediate, and large Ta concentrations 
the ground state is c(2 x 2)-AFM, p(2 x 1)-AFM, and FM, 
respectively (Fig. [5]). The qualitative agreement between 
FLAPW-VCA and TB-LMTO-CPA proves that the VCA 
is a good approximation to the CPA for the treatment of 
the alloy. Interestingly, for intermediate compositions the 
DLM and the p(2 x 1)-AFM states are degenerate, within 
the computational accuracy. Note that the DLM state is 



often encountered as a spin-glass-like ground state in dis- 
ordered bulk alloys, such as fee Ni .8oMn .2o [HI- All im- 
portant conditions for the spin-glass (SG) arrangement, 
namely competing ferro- and antiferromagnetic exchange 
interactions accompanied by chemical and/or topological 
disorder [23| , are fulfilled in the present case owing to the 
random Ta-W substrate. At around x — 85% Ta concen- 
tration, the SG, FM, and p(2 x 1)-AFM states are de- 
generate. This offers the great perspective of imprinting 
magnetic information as nanoscale dots. E.g., cooling 
down the dot, after local heating at elevated tempera- 
tures, with or without an external magnetic field allows 
to systematically select the FM or the zero-magnetization 
SG state. 

In summary, we introduced here a new model system 
to study magnetic phase transitions in low-dimensional 
magnets allowing to sweep through a series of magnetic 
phases without changing the structure or chemical com- 
position of the magnet. The system offers the possibility 
to store information in terms of FM dots in an otherwise 
zero-magnetization state matrix. 
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